Mean-Field Analysis of Spinor Bosons in Optical Superlattices 
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We study the ground-state phase diagram of spinless and spin-1 bosons in optical superlattices using a Bose- 
Hubbard Hamiltonian that includes spin-dependent interactions. We decouple the unit cells of the superlattice 
via a mean-field approach and take into account the dynamics within the unit cell exactly. The system supports 
Mott-insulating as well as superfluid phases. The transitions between these phases are second-order for spinless 
bosons and second- or first-order for spin-1 bosons. Anti-ferromagnetic interactions energetically penalize high- 
spin configurations and elongate all Mott lobes, especially the ones corresponding to an even atom number on 
each lattice site. We find that the quadratic Zeeman effect lifts the degeneracy between different polar superfluid 
phases leading to additional metastable phases and first-order phase transitions. Finally, we show that an energy 
offset between the two sites of the unit cell induces a staircase of single-atom tunneling resonances which 
surprisingly survives well into the superfluid regime. 

PACS numbers: 03.75.Lm, 67.85.Fg, ()3.75.Mn, ()5.3().Rt 



C/3 

a- 



o 



0^■ 
O 



X: 



I. INTRODUCTION 

Ultracold atom experiments offer the unique opportunity to 
study quantum many-body effects in an extremely clean and 
well-controlled environment. In contrast to most condensed 
matter systems they are characterized by the absence of disor- 
der and other imperfections and are highly controllable. This 
is why they have been proposed as quantum simulators Ol-dl] . 

One of the most prominent achievements in this direction 
has been the observation of the quantum phase transition be- 
tween a Mott insulating and a superfluid phases of ultracold 
atoms in an optical lattice (^-^. In this situation, the ratio of 
the tunneling strength between the lattice sites and the on-site 
interaction determines the qualitative behavior of a quantum 
gas. If the on-site interaction dominates and the filling is com- 
mensurate the quantum gas is in the Mott-insulating phase. 
On the contrary, if the tunneling amplitude is sufficiently large 
the system becomes superfluid. 

For atoms trapped in a magneto-optical trap the spin degree 
of freedom is frozen and the atoms become effectively spin- 
less. If, however, the quantum gas is trapped by purely opti- 
cal means, the atoms keep their spin degree of freedom and 
the order parameter describing the superfluid phase becomes 
a spinor. These systems were first studied in Refs. |0, [Ml and 
recently reviewed in Ref . . 

Due to their spin-dependent interactions ultracold spinor 
quantum gases in optical lattices offer the possibility to model 
mesoscopic magnetism [26]. They are well described by the 
Bose-Hubbard model, but the spin-dependent effects alter the 
system in a qualitative and quantitative way lflol - [T5ll . The 
phase boundaries between superfluid and a Mott-insulating 
phases are shifted, and for certain atom-configurations the 
phase transition is no longer second- but first-order. This is 
a consequence of the additional spin-dependent on-site inter- 
action. If this interaction is anti-feiTomagnetic, atomic sin- 
glets are energetically favored and the Mott-insulating phase 
is stable in some parameter ranges where the system is super- 
fluid for spinless atoms. The occuiTence of first-order phase 
transitions enables metastable phases and the system shows 
hysteretic behavior. 




FIG. 1 : (Color online) Potential landscape of an optical superlattice 
and parameters of the Bose-Hubbard model {2} for spinless ultracold 
atoms in optical superlattices. The dots depict the atoms in the super- 
lattice potential V{x). e is the energy offset between the two sides of 
the double well, ti (te) is the intra- (inter-)well tunneling amplitude, 
and U is the strength of the on-site interaction. 



In this paper we investigate ultracold bosons loaded into 
optical superlattices. These systems have already been stud- 
ied theoretically lfl6[ [T7tl and experimentally flsl - Bli] . We take 
into account the dynamics in the double wells exactly and in- 
clude the tunneling between neighboring unit cells via a mean- 
field ansatz. The system supports Mott-insulating phases as 
well as superfluid phases Iq3l l23ll . The former are charac- 
terized by a fixed number of atoms per unit cell. For spin-1 
atoms we find that some of the phase transitions become first 
order similar to the case of usual period- 1 lattices lITill . 

We include the effects of magnetic fields by using an ef- 
fective Hamiltonian which includes a quadratic Zeeman shift. 
For anti-ferromagnetic interactions magnetic fields break the 
degeneracy between different polar superfluid phases. This 
leads to new classes of metastable phases and changes 
the phase boundaries significantly. In the ferromagnetic 
case magnetic fields cause first-order phase transitions and 
metastable phases. These results apply to spin-1 atoms in su- 
perlattices as well as in usual lattices. 

Finally, we examine single-atom tunneling resonances in 
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FIG. 2: (Color online) Sketch of the two-dimensional period-2 su- 
perlattice. The filled circles depict the left wells of each unit cell, 
the open ones the right wells. The tunneling amplitude within the 
unit cells is ti. There are three different tunneling processes between 



neighboring unit cells, tl}\ t^^\ 



and t 



(3) 



superlattices. These are known from isolated double-well po- 
tentials ll24l - t26ll and are generalized to superlattices in this 
work. 

The remainder of this paper is organized as follows. In 
Sec. |II]we introduce the mean-field Hamiltonian for spinless 
bosons in optical superlattices. We discuss methods to treat 
this Hamiltonian and present the phase diagram in Sec. Ill Al 
Section HTB] treats the phenomenon of single-atom resonances 
for spinless atoms in optical superlattices. In Sec.|lII]we gen- 
eralize the Bose-Hubbard Hamiltonian of Sec.|II]by including 
spin-dependent interactions. In Sec. IIII Al the phase diagram 
of spin-1 atoms is examined. We discuss the novel aspects 
due to the spinor nature of the order parameter and point out 
differences in the single-atom resonances. Finally, we include 
magnetic fields in Sec. IIII B I which enhance spin-dependent 
effects and lead to additional metastable phases. 



n. SPINLESS BOSONS 

Optical lattices are arrays of weakly connected optical 
micro-traps created by interfering counter-propagating laser 
beams. When the lattice is superimposed with a second stand- 
ing laser beam with half the wavelength an optical superlattice 
is created (see Fig.[T]). The unit cell of a period-2 superlattice 
is a double-well potential lll"8l - (2lll . 

Ultracold bosons in sufficiently deep optical lattices can be 
described by the Bose-Hubbard Hamiltonian |5, 27]. Here we 
examine the physics of bosonic atoms in optical superlattices 
when the overall atom density is chosen such that there are 
a few atoms per double well. If one neglects tunneling be- 
tween neighboring unit cells the atoms in each unit cell can be 
described by the Hamiltonian 



U 
2" 



k=L.,R 

+ e{nL- ur) - y.{nL+ fin) , 



where L (U) and R {u') are bosonic annihilation (creation) 
operators for atoms in the left or right well, {fiji) is the 
atom number operator at the left (right) site. U is the on-site 
interaction and ti is the tunneling strength between the sites of 
the double well. The energy offset between the sites is given 
by e and the chemical potential is fi (see Fig. [TJ. The pa- 
rameters can be tuned by changing the intensity and the phase 
difference between the counter-propagating laser beams; it is 
possible to tune the system from the regime of strong tunnel- 
ing {ti ^ U) to the regime of weak tunneling (ti ^ U). 

The Hamiltonian of an array of connected double-well po- 
tentials includes tunneling between neighboring unit cells. We 
choose the configuration as shown in Fig.|2]where there are in 
general three different inter- well tunneling amplitudes, ti^\ 

(2) (3) 

te and te ■ It turns out that the our results depend only 
weakly on the differences among the inter-well tunneling am- 
plitudes. This is why we will assume te = ii^^ = i?'' — t''J^\ 
Within the mean-field approximation which has been de- 
veloped in Refs. ll22ll28ll . our Hamiltonian reads 

H = Ho-te {4>rL'^ + 0L^t _^ 2z0fl.7?t + 2z(j)LL^ 

- 4>r4>*l - z(t)R(t)*R - z(t)l(j)L + h.c)j (2) 

where we introduced the mean-field parameters (t)n = (R), 
(f)L = (i), and 2 = 1 for 2D lattices and z — 2 for 3D lattices. 
The Hamiltonian treats the internal degrees of freedom of each 
unit cell exactly and approximates the tunneling between the 
unit cells via a mean-field ansatz. This approach is expected 
to give satisfactory results if the tunneling inside the unit cells 
is stronger than the tunneling between the unit cells (i.e., te < 
ti), otherwise correlations between neighboring double wells 
would be stronger than correlations within the double wells 
and should not be neglected. 

The system is in the Mott-insulating phase if (j)L — (f>R — 
and in the superfluid phase if (f>L ^ ^ (f>R. In the latter 
case the number of superfluid atoms and on the left 

and right site is given hy $ = {4)l,4'r) = {\J~^, \f^)- 

There are several possibilities to treat the Hamiltonian (|2|l. 
For a given set of parameters {pi,e,ti,te, U} the task is to 
find the self-consistent values of — i r 

/,(o)i£.|„/,(o)\ 



</)L = {,p'^"\L\i;p, and 

denotes the ground state of 

the Hamiltonian for a given order parameter (p. On the one 
hand the self-consistent values are fixed points of the map 



bR^{,p'^>\R\i^'i'>) where 



(3) 



where the index i refers to the ith step in the iterative pro- 
cedure used to find the self-consistent value of the order pa- 
rameter. On the other hand the self-consistent values of (p 
correspond to the local extrema of the energy functional 



/,(o)\ 



(4) 



(1) 



and its local minima correspond to stable fixed points of the 
map (O which can be found by the iterative procedure. 
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FIG. 3: (Color online) Ground-state energy E[(f>L , (pR,] of the Hamil- 
tonian ([2} as a function of the order parameter — {(j>L, 4'r\- For 
a symmetric unit cell e = 0, we have <j)L = (i>R ~ (j>- The red 
line _B[(/)]'^°" corresponds to the Mott-insulating phase {n/U = 0.2, 
ti/U = 0.05, and te/U = 0.005), the blue line E[<f)Y^ to the super- 
fluid phase ip/U = 0.5, U/U = 0.22, and te/U = 0.022). 
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FIG. 4: Phase diagram of spinless bosons in a superlattice described 
by the Hamiltonian In the atomic limit ti = = G the Hamilto- 
nian is diagonal in the Fock basis and supports only Mott-insulating 
phases. The blue lines mark the phase boundaries and (ni,, n^) de- 
notes the occupation of the left and the right site in the double well. 




FIG. 5: (Color online) Phase diagram of spinless bosons in a two- 
dimensional superlattice described by the Hamiltonian l|2ll. We plot 
the critical internal tunneling amplitude ti as a function of chemical 
potential and energy offsets fort, = IQte. In Fig. g] we show a 
cut through this 3D plot at = te = and in Fig.|6]at ti/U = 0.05. 
The edge at e = of the phase diagram reveals the contraction of 
Mott lobes to loops at integer values of ^/U . 
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FIG. 6: (Color online) Phase diagram of spinless bosons in a 
two-dimensional superlattice described by the Hamiltonian ^ for 
ti/U = lOte/f/ = 0.05. The blue areas depict the superfluid phase. 
The dashed line corresponds to the parameters chosen for Fig. [7] 



For the Hamiltonian (|2]l there are only two classes of energy 
functionals (see Fig.O: those with only one local extremum 
at (f)L ^ 4iR = Q corresponding to a Mott-insulating phase, 
and those with a second extremum at 7^ 7^ which is 
the global minimum and corresponding to a superfluid phase. 

This enables us to distinguish the Mott and superfluid quan- 
tum phases with minimal numerical effort. We only have to 



calculate the ground-state energy for = and in its proxim- 
ity (/) ss 0. If E[<l) « 0] — E[(t) = 0] is positive, the system is 
Mott-insulating; if it is negative the system is superfluid. 
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A. The phase diagram for spinless bosons 



In this section we determine the ground-state phase diagram 
of the Hamihonian (|2|i. For a chemical potential /i, an energy 
offset e, and a given ratio of the tunneling amplitudes te/ti 
we calculate the critical tunneling amplitude above which 
the system is superfluid. In the following we use the on-site 
interaction U as the unit of energy. 

In Fig.|4]we plot the ground state as a function of chemical 
potential /i and offset e in the atomic limit ti = = 0. In this 
case, the Hamiltonian ^ is diagonal in the Pock basis and the 
system supports only Mott phases {riL^nii) characterized by 
the number of atoms in the left nj^ and right well n^. Each 
of the diamonds in Fig.|4]corresponds to one Fock state, i.e., a 
fixed particle number in the unit cell as well as a fixed particle 
number in the left and the right site of each unit cell. When we 
increase /i/C/ for fixed energy offset e/U the number of parti- 
cles in the unit cells increases while the ratio between left and 
right remains similar When we increase e/U for fixed chem- 
ical potential n/U the atom number is constant but the atom 
distribution within the unit cells changes. Because we set the 
tunneling to zero this happens non-continuously. Figure |4] is 
mirror-symmetric along the e = axis, i.e., when e — > — e the 
atom number distribution is inverted, (n^, nji) [na, ul). 

In Fig.|5]we plot the critical tunneling strength ti at which 
the system becomes superfluid as a function of the chemical 
potential ji and the offset e for ti/t^ — 10. It is convenient 
to pick a fixed ratio of ti / in order to obey the constraint 
te < ti. For fixed energy offset e we recover Mott lobes, 
which are familiar from the case of a usual lattice |4], although 
the Mott phase for atoms in superlattices is characterized by a 
fixed integer atom number per unit cell, i.e. n — {n^ + fiR) 
where n is an integer number When e/U has an integer value 
the Mott lobes corresponding to an odd atom number per unit 
cell contract to Mott loops and if e/U has an half-integer value 
the lobes corresponding to an even atom number contract to 
loops [22]. As the energy offset e is varied, the size of the Mott 
lobes changes and they constitute tubes of fixed integer atom 
number per unit cell. The base of the plot (i.e., the ti ^ t^ = 
plane) shows the diamond structure given in Fig.|4] The nodes 
of the diamonds are special: these are the values of the energy 
offset e where the lobes contract to loops, i.e. the Mott tubes 
touch the ti = <e = plane only at one point. 

Figure |6] presents a cut through Fig. |5] at ti/U = 0.05, 
showing the Mott insulating phases in white and the superfluid 
phases in blue. The Mott diamonds of Fig. |4]are connected 
for non- vanishing tunneling amplitudes. This means that the 
quantum numbers {11^,71 ji) change continuously when e is 
varied and the Mott insulating phases are characterized by one 
number n = ul + nu, the total number of particles per unit 
cell. In Fig. |6]the chosen tunneling amplitudes are too large 
to see the connections between the Mott diamonds for n > 5, 
nevertheless the quantum numbers (n^, n^) are not fixed to 
integer values for these Mott phases either. 
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FIG. 7: (Color online) Mean particle number in the right well (nn) 
(solid) and standard deviation Ann ~ \/ {'"■%} ^ ('^-r)^ (dashed) as 
a function of the energy offset e for Hamiltonian (O with jj./U — 1 
andti/U = lOte/U — 0.05. The inset shows the mean total particle 
number in the double well (n) = ifiL +nR). The standard deviation 
is multiplied by 5. 
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FIG. 8: (Color online) Mean particle number in the right well (fin) 
(solid) and standard deviation Ann ~ \/ ('^/j) ^ ('^-r)^ (dashed) as 
a function of the energy offset e for Hamiltonian ^ with n/U = 
1.2 and ti/U — lOte/U = 0.1. The inset shows the total particle 
number in the double well (n) = {fiL +nR). The standard deviation 
is multiplied by 5. 



B. Single-atom resonances within superlattices 

If the offset e is changed with the other parameters /i, ti, 
and fixed, the atom number distribution within the double 
wells becomes asymmetric. Due to the finite on-site interac- 
tion the atom numbers do not change proportionally to e but in 
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FIG. 9: Fraction of condensed particles in the right well n''^/n^^ ([5j 
for spinless atoms described by the Hamiltonian The parameters 
are chosen such that the system is superfluid over the full range of 
e/U £ [—1.8, 1.8]. The chemical potential for all curves is n/U — 
1.25. The solid line corresponds to ti/U = lOte/U = 0.15, the 
dashed line to ti/U — lOte/U — 1 and the dashed-dotted line to 
ti = te — <x. The dotted line shows the atom number in the right 
well as a function of efoi ti/U = 10te/U = 0.15. 



steps, which are called single-atom resonances. For isolated 
double-well potentials these bosonic staircases were theoreti- 
cally predicted [24, 26, 29] and experimentally detected jzsll . 

Within the Mott-insulating phase the superlattice decom- 
poses into an array of isolated double wells. Thus, it is possi- 
ble in this regime to observe the same bosonic staircases as in 
the case of single double-well potentials. 

Figure [T] shows the mean atom number in the right well 

\2 along the 



(77, fl) and the standard deviation {nj^j — [na) 
dashed line shown in Fig. |6] The system is Mott-insulating 
for a large parameter range (ii/U = 1, ti/U = 0.05, te/U — 
0.005, and -1.9 < e/U < 1.9). The inset shows the total 
number of atoms in the double well (h) = {ul + fin) which 
changes discontinuously aie/U « ±2 signaling that the sys- 
tem leaves the Mott-insulating phase. This causes a discon- 
tinuity in the function {hn). The steps are equidistant since 
the difference in the on-site interaction between one and two 
atoms and between two and three atoms (and for higher atom 
numbers) are the same, namely equal to the on-site interaction 
U . The spacing between two steps is Ae — U, i.e., the steps 
occur when the energy offset is large enough to compensate 
the on-site interaction. 

In Fig.|8]the expectation value {hn) is plotted along a curve 
with stronger tunneling amplitudes and a shifted chemical po- 
tential as compared to Fig. [T] Along this curve the system is 
mostly superfluid. This can be seen from the inset of Fig. |8] 
For —2U < e < 2U the particle number per double well {n) 
is not constant anymore. Although the atom staircase looks 
similar to the one in Fig. [T] there are significant differences. 
The fluctuations in the particle number An^ are greatly in- 
creased and additional maxima appear. These number fluctu- 



ations can be measured in experiments fz?) . 

Single-atom resonances can also be seen in the density of 
condensed atoms. The number of condensed atoms per site is 
connected to the value of the order parameter via the relations, 



\{L)\- 



and n^? 



I (i?) I . The ratio of condensed atoms 



in the right well to the number of condensed atoms in both 
wells. 



1(^)1^ 



n^* \{LW + \{R)\ 



2 ' 



(5) 



is plotted in Fig.|9]along several paths in parameter space. The 
solid line corresponds to ti — lOtg — 0.15U and shows step- 
like behavior. The dotted line shows the total atom number in 
the right well {hn) for the same tunneling amplitudes. Note 
that the steplike structure is more pronounced for the density 
of condensed atoms than for the total atom number. For higher 
tunneling rates the staircase structure disappears (see dashed 
line with ti — lOt^ = lU). In the limit of infinite tunneling 
amplitudes the ratio n^/j^/n'^^ does not depend on e because 
there is an infinite amount of atoms in the unit cell and chang- 
ing e/U hy one moves only one atom from one site to the 
other. This analysis of the asymmetry of the superfluid den- 
sity also helps to understand the connection between 0^ and 
(f)ii. The ratio of the two order parameters (p^ and (f>ii does not 
only depend on the energy offset e but also on the tunneling 
amplitudes, the chemical potential, and the on-site interaction. 



ni. SPINOR BOSONS 

The Bose-Hubbard model of Eq. ^ can be generalized to 
spin-1 atoms by taking a spin-dependent on-site interaction 
into account 1 10]. To obtain a double-well mean-field Hamil- 
tonian for spin-1 atoms in optical superlattices we have to re- 
place in Eq. (|2]i the annihilation and creation operators, as well 
as the order parameters, by 3-component vectors according to 
the three hyperfine projections for spin-1 bosons and add a 
term containing the spin-dependent on-site interaction. The 
resulting Hamiltonian is 

H =Y ^kifik - 1) - U (l^ R + /i.e.) 

k=L,B. 

+ e{fiL- flu) - fJ.{hL+ flu) 

U. ^ 



k=L,R 



ie • lV (/-L • R' 

$R- $L- z4>L ■ <P*L - Z(f)*ji ■ (t)R + h.C. (6) 



whereL' — {lJ, ij, L^-^} is a vector containing the creation 
operators for the left well, i.e., L„j creates an atom in the mth 
hyperfine state in the left well. Similarly, L consists of an- 
nihilation operators of the left well and R (R^) are annihila- 
tion (creation) operators for the right well, — J2a ^cr^cr 
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FIG. 10: (Color online) Ground-state energy E[(f>L,4>R] of the 
Hamiltonian l[6]l as a function of the order parameter (j). The unit 
cell is chosen to be symmetric e = 0, i.e., we have (f)L = (fiR ~ (j), 
we have chosen (po — 0, and 01 — (f)-\ due to symmetry constraints. 
The blue line corresponds a point in parameter space 

{H/Uo = 0.25, U/Uo = 0.35, U/Uo = 0.035, U2/UQ = 0.04 , and 
£ = 0) where there is a metastable Mott phase in addition to the sta- 
ble superfluid phase. The green line E[<f)\'^'^^ corresponds a point 
within the metastable superfluid phase ip/Uo = 0.25, ti/Uo — 0.3, 



te/Uo = 0.03, U2/U0 = 0.04, and e = 0). E[, 
0.09(/o to show the two curves in the same plot. 



is shifted by 



hff = X)cr R^jRcrj is the atom number at the left (right) 
site. The annihilation and creation operators obey the canon- 
ical commutation relations [Li,Lj] = [Ri,R}^] — and 
ti _ r?. £.ti 



[U,R]] = 0. Si = Haa'Ll'iaa'L^' iS the 

total spin on the left site and the total spin on the right site 

is Sij — X^ctct' R-t'icrcr'Ra', wherc Taa' are the usual spin- 
1 matrices. The dimensionality of the array is contained 
in the parameter z; for 2D lattices z — 1 and for 3D lat- 
tices z = 2. The vectors (fi, = {'/'l \ '^l'''' "^l 
$R = {'p'^R 1 'i'^R J } contain the six order parameters of 
the Hamiltonian (|6]l. Note that the system is rotationally sym- 
metric and 01 = 0_i for both 0^ and $ft. 101 ■ 

The term proportional to U2 describes spin-dependent con- 
tact interactions: in the case of anti-ferromagnetic interactions 
U2 > (e.g. ^■^Na) it penalizes non-zero spin configurations 
while it favors high-spin configurations in the case of ferro- 
magnetic interactions U2 < (e.g. *^Rb). Whereas the ratio 
t/Uo can be controlled with the intensity of the laser beams, 
the ratio U2/UQ depends on the spin-2 and spin-0 scattering 
lengths of the spin-1 atoms [7], e.g. U2/U0 = 0.04 for ^^Na. 

The Hamiltonian ^ has a much richer phase diagram than 
the Hamiltonian (|2]i. In addition to Mott-insulating and super- 
fluid quantum phases, the spin-1 Bose-Hubbard model gives 
rise to metastable quantum phases. This can be seen by look- 
ing at the energy functionals E[(j)L, ipn] of the Hamiltonian 
(|6]l. In addition to those shown in Fig.|3] two other classes of 



FIG. 11: (Color online) Energy functional E[(j>L,4'R] for the 
Hamiltonian ^ with anti-ferromagnetic interactions (ti/Uo = 0.3, 
te/Uo = 0.03, U2/U0 = 0.04, fi/Uo = 0.25, e = 0). Because the 
unit cell is chosen to be symmetric the order parameters of the left 
and the right are equal {(j) — cj>L = 4'r) and the corresponding in- 
dices are suppressed. Because of the rotational symmetry <jf>i = <j)--i. 



energy functionals arise for anti-ferromagnetic spin interac- 
tions U2 > 0, see Fig.[TOl These energy functionals have two 
local minima and an iterative procedure similar to the one de- 
scribed in SeclUdoes not lead to an unique value of the order 
parameter but depends on the starting point. When the start- 
ing point of the iterative procedure is chosen close to zero, one 
finds the minimum at 0l = = 0. If one starts at a value 
beyond the maximum separating the two minima one obtains 
the second minimum corresponding to a superfluid phase. The 
global minimum of the energy functional determines the sta- 
ble quantum phase of the system. The other one corresponds 
to a metastable phase. 

Thus, energy functionals such as the ones in Fig. [TOl signal 
metastable phases, first-order phase transitions and hysteric 
behavior of the system; they do not allow the same analysis 
as the spinless case. The stability analysis of the (f> — fixed 
point does not answer the question, if there is a second stable 
fixed point and if it is energetically lower or higher. To deter- 
mine the quantum phase we numerically calculate the energy 
functional and analyze its local minima. 

Due to the spinor nature of the order parameter additional 
properties of the superfluid phases arise. The spin-dependent 
interaction changes the symmetry of the energy functional 
E[(Pl,(I}r] in the (j)o-4>i plane for 0^ as well as For 
1/2 = the energy functional is rotationally symmetric around 
<l>i = <t>o = 0. For anti-ferromagnetic interactions U2 > the 
bosons form a polar superfluid, i.e. the spin-dependent inter- 
action energy is minimized by (Sl) = {^r) = 0j|7[l- There 
are two different classes of polar order parameters fTstl, one is 
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the transverse polar state 




ci 



the other one the longitudinal polar state 



= Co 



(7) 




(8) 



where cq and ci are real numbers. Here we omitted the in- 
dex labeling the left and right site because the structure of the 
order parameter is the same for both of them. 

In Fig. [TT]we plot the energy functional E[(I)l,(Pl] as a 
function of the order parameter for a symmetric unit cell and 
a point in parameter space at which there are two local min- 
ima corresponding to a longitudinal and a transverse polar su- 
perfluid phase as well as a local minimum signaling a sta- 
ble Mott insulating phase. Since the unit cell is chosen to be 
symmetric (e — 0) the order parameters in the left and the 
right site are the same (0^ = (pfj). The two minima at a 
non-vanishing order parameter are degenerate, because both 
correspond to (Sl) = (Su) = and therefore suffer the same 
spin-dependent energy shift of the on-site interaction. Due to 
the special form of the two superfluid phases given in Eq. Q 
and Eq. (|8]l the superfluid minima are always on the ipi = 
and the 0o = axes, respectively. This justifies why we chose 
(j)o = in Fig. [TO] The additional minimum in Fig. [TT] at 
(j)o = (pi — (p-i ~ is the global minimum and corresponds 
to the Mott-insulating phase; the two degenerate minima cor- 
responding to a non-vanishing order parameter belong there- 
fore to two degenerate metastable phases. 

In the ferromagnetic case (U2 < 0) there is only one su- 
perfluid order parameter which is given by ((/>i, (f)o,(f>^i)^ = 
c(l,\/2,l)^, where c is a real number Again, we suppressed 
the index labeling the left and right site because it is the same 
for both. In the ferromagnetic case, the spin-dependent inter- 
action in Eq. (|6]l has the same sign as the tunneling term and 
therefore does not create metastable quantum phases. 



A. The phase diagram 




FIG. 12: (Color online) Phase boundaries of spin-1 bosons in a two- 
dimensional superlattice for ti = = 0, described by the Hamil- 
tonian {Sjl. The red lines denote deviations in the phase boundaries 
relative to the spinless case (see Fig. |4](, black solid lines to phase 
boundaries which are not changed and black dotted lines to shifted 
phase boundaries of the spinless case. Each Mott diamond is labeled 
by its atom number configuration (n_L,n_R) (in black) and the en- 
ergy penalty due to spin-dependent interactions (in red below). The 
green diamonds correspond to odd-even particle number configura- 
tions, the yellow ones to odd-odd, and the blue ones to even-even 
configurations. 
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TABLE I: Diagonal elements of the Hamiltonian (|6} for ti — te — 
in Fock space (i.e., \nL,nR)). Because we are interested in ground- 
state properties for anti-ferromagnetic interactions we choose for 
each atom number configuration the smallest spin configuration. For 
ferromagnetic interactions the highest spin configuration is energeti- 
cally favorable and the table is thus changed. 



In this section we calculate the phase diagram of spin-1 
atoms in superlattices. We focus on the differences to the spin- 
less case which was discussed in Sec. Ill Al 

For ti ~ t,, — the Hamiltonian ^ is diagonal in the 
Fock basis and the system supports only Mott phases. Simi- 
lar to the spinless case the Mott phases are characterized by 
{nL,nji) and the boundaries between Mott phases carrying 
a small number of atoms can be calculated directly from the 
collection of eigenstates and eigenenergies in Table U 

In Fig. [12] we show how the phase boundaries are shifted 
compared to the spinless case for anti-ferromagnetic interac- 
tions. Again, there is a diamond-shaped structure as in Fig.|4] 
The Mott diamonds of Fig. |4] increase or shrink depending 



on their atom number configuration. For anti-ferromagnetic 
spin interactions, the strength of the on-site interaction de- 
pends upon the parity of the atom number at each lattice site. 
An even number of spin-1 atoms allows the formation of a 
spin singlet, i.e. vanishing total spin per site, which minimizes 
the on-site repulsion. Odd atom numbers are penalized, be- 
cause the spin-singlet wavefunction is anti-symmetric for an 
odd atom number and thereby ruled out by symmetry con- 
straints. Diamonds corresponding to an even particle number 
in the left as well as the right well are favored and diamonds 
corresponding to an odd-odd configuration are penalized. 
The boundary between two Mott diamonds is shifted only if 



8 




FIG. 13: (Color online) Phase diagram for spin-1 atoms in optical 
supeiiattices (U2/U0 ~ 0.04, e = 0, and ti — lOte). The shaded re- 
gions denote Mott-insulating phases. The dashed lines are the phase 
boundaries for metastable phases and the dotted lines are the phase 
boundaries for the spinless case U2 = 0. The regions in which a 
metastable Mott phase coexists beside the superfluid [SF] phase is 
marked with MM; MSF denotes regions where metastable superfluid 
phases exist alongside the Mott [Mott] phase. The Mott lobes are 
labeled according to the total atom number n per double well. 




1 2 s/Uo 



FIG. 14: (Color online) Bosonic staircase for spin-1 atoms described 
by the Hamiltonian ^ along the curve /i/f/o — 1.4, ti/Uo ~ 0.05, 
te/Uo = 0.005, U2/U0 = 0.1 and £/[/(, G [-2.5,2.5]. The thick 
solid line shows the particle number in the right well (hr) and the 
dashed line the standard deviation Ann = \/ {n^} — {tir)^. The 

thin solid line shows Sr, the quantum number of the square of the 

" 2 

spin in the right well, i.e. (Sr) = Sr{Sr + 1). The total particle 
number {ul + ur) = 4 over the full shown range of e. The standard 
deviation is multiplied by 4. 



the spin-dependent energy penalty (or bonus) is different for 
them (compare Table The value of this shift depends lin- 
early on [/2. In the anti-ferromagnetic case, a phase boundary 
is either unshifted or shifted by a constant amount. This is be- 
cause at each phase boundary the atom configuration changes 
only by one atom and the only possible ground-state spin con- 

figurations at each lattice site are spin singlets ((S ) = 0) 

^ 2 

and total spin equals one ((S ) = 2). This is the reason 
why the straight lines in Fig. |4]are preserved in case of anti- 
ferromagnetic interactions. Thus, to determine the shift of the 
phase boundaries it is enough to examine an example, say 
the phase boundary between the Mott diamonds containing 
one atom in the left well and the ones containing two atoms 
in the left well. The phase boundary follows for positive e 
the path along e ~ — fi + Uq — 2U2, which can be seen by 
setting Ell — E12 or Eiq = i?2o- The shift of the phase 
boundary is therefore A ~ \/2U2- The diamonds correspond- 
ing to an odd number of atoms in the double well (green di- 
amonds in Fig. \T2\ change their size from 1/2 (in units of 
Uq) to 1/2 — A'^/Uq. Double wells caiTying an even num- 
ber of atoms allow even-even configurations (blue diamonds 
in Fig.[T2]with area (l/\/2 + A/C/q)^) and odd-odd configura- 
tions (yellow diamonds in Fig . [12] with area ( 1 / — A/Uo)^)- 
The full ground-state phase diagram for anti-ferromagnetic 
interactions, a symmetric double well and ti = lOig is shown 
in Fig. [T3] We choose U2/U0 = 0.04 corresponding to 
^■^Na |3 1] . Spin-dependent interactions lead to elongated IMott 
lobes. In general, the spin configuration is higher in the su- 



perfluid phase than in the IVIott phase and this leads to an en- 
ergy penalty (see Hamiltonian (|6]l). Whenever there is an even 
number of atoms in a lattice site this effect is strongest, be- 
cause an even number of spin-1 bosons can form spin singlets 
(see above). In Fig.[T3]the Mott lobe containing four atoms is 
therefore significantly enlarged, since this Mott phase corre- 
sponds to two atoms in the left as well as the right site (the unit 
cells are chosen to be symmetric, i.e., e — 0). The Mott lobe 
n = 2 is significantly enlarged for e/Uo = ±1 (not contained 
in Fig.[T3]l. because the atoms pair up on the left (or right) site 
and form spin singlets. The phase transitions between Mott 
lobes corresponding to an even number of atoms in the dou- 
ble well are of second order whereas the others are first-order 
phase transitions. For smaller values of U2/U0 all phase tran- 
sitions become first order at a tricritical point, in contrast to 
spin-1 atoms in usual lattices for which the boundary of the 
Mott phase with one atom per site is always a second-order 
phase transition 1 14]. 

To illustrate the phase diagram and to compare it to the case 
of spinless bosons we include a bosonic staircase for spin- 
1 atoms with anti-ferromagnetic interactions in Fig. [141 In 
this figure the occupation number of the right site of each unit 
cell {hji) is plotted as a function of the energy offset e. We 
choose parameters so that the system is Mott-insulating for 
—2.5 < e/Uo < 2.5. The comparison with Fig.|2]shows that 
anti-feiTomagnetic interactions shift the steps and make them 
unequally wide. The step corresponding to two atoms in the 
left site and two atoms in the right site allows the formation 

of spin singlets ((S"l) = (S^;) = 0) in both sites. This is 
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FIG. 1 5 : (Color online) Energy functional E[<I)l,(I>r] for a symmet- 
ric unit cell (e = 0) and anti-ferromagnetic interactions U2/U0 ~ 
0.04 in the presence of a magnetic field q/Uo = 0.002. The parame- 
ters are chosen identical to Fig.[TT](ti/i7o ~ 0.3, te/Uo = 0.03 and 
fi/Uo = 0.25). Because the unit cell is chosen to be symmetric the 
order parameters of the left and the right site are equal (<I>l ~ <I>r) 
and the corresponding indexes are suppressed. Because of the rota- 
tional symmetry = 0-1. 



energetically favorable compared to the case of an odd number 

of atoms on both sites, which form states with (S^) = (S^) = 
2 (this is analogous to the changed size of the diamonds in 
Fig.IIl. 

B. Magnetic fields 

Finally, we want to examine the effect of weak magnetic 
fields on the phase diagram of spin-1 atoms in optical lattices. 
Spinor Bose-Einstein condensates are realized in optical traps, 
because strong magnetic fields would break the degeneracy of 
the different hyperfine levels and align the atom spins, thus 
creating a condensate of effectively spinless bosons. Never- 
theless it is interesting to include weak magnetic fields in the 
calculations because magnetic shielding is needed but can be 
done only to a certain degree. To incorporate the effect of a 
magnetic field in the Bose-Hubbard model for spin-1 bosons 
one can add a term to the Hamiltonian given in Eq. (|6]l which 
includes the Zeeman shift of the energy levels |30]. The first 
contribution of a magnetic field B = (0, 0, B) chosen along 
the 2-axis (which is our quantization axis) is a linear Zeeman 
shift, 

i=L,R a 

where p — qhbB and fiicr is the particle number operator 
for the ith site that gives the number of bosons in the mth 
hyperfine state. Although the linear Zeeman shift alters the 
energy levels considerably, it can in general be neglected. The 
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FIG. 16: (Color online) Phase diagram for spin-1 atoms in opti- 
cal superlattices with a magnetic field described by the Hamiltonian 
d {U2/U0 = 0.04, e = 0, ti = lOte, and q/Uo = 0.02). 
The shaded regions denote Mott-insulating phases. The dashed 
lines are the phase boundaries for metastable phases. There are 
no metastable quantum phases in the Mott-insulating region, but 
additional metastable phases in the superfluid one. The quantum 
phase (a) denotes a superfluid region where there are two metastable 
phases. There is a metastable superfluid phase of the transverse po- 
lar state and energetically higher a metastable Mott phase. In the 
quantum phase (b) there are also both metastable phases but the 
metastable Mott phase of the transverse polar state is energetically 
lower than the metastable superfluid phase. The quantum phase 
(c) contains only one metastable phase of the transverse polar state 
which is a Mott phase. 



interactions within the atomic cloud do not change the overall 
magnetization [s^, [Ht] and the linear Zeeman shift can be 
gauged away. For a given system with fixed magnetization 
the main impact of a magnetic field is given by the quadratic 
Zeeman effect, 

Hqz = 9 X! X!"^^'^^*'^' ^^^^ 

i=L.,R (J 

which is added to the Hamiltonian (|6]l. The magnitude of the 
quadratic Zeeman shift is given by (7 = QqB^, where e.g. qo = 
h X 390 Hz/G^ for sodium HI. 

The quadratic Zeeman effect affects the phase diagram 
considerably. The local minima of the energy functional 
E[(t>Li <t>R\ belonging to transverse and longitudinal polar su- 
perfluids are no longer degenerate (see Fig. [T5] l. For pos- 
itive q the longitudinal superfluid states are always ener- 
getically favored, for negative q the transverse ones. For 
anti-ferromagnetic spin interaction new classes of metastable 
quantum phases arise, see Fig. [16] It is important to notice, 
that even very weak magnetic fields (q/Uo — 0.002 in Fig.fTSl) 
change the properties of the ground state energy functional 
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substantially. In Fig. [16] we choose the same parameters as 
in Fig. [TTIjust adding a very weak magnetic field. The mag- 
netic field causes a quantum phase transition from the super- 
fluid phase to a Mott insulating one since the minimum at 
(j)i — = Q, (t)Q — 0.55 is now the global one. Addition- 
ally, there are now two metastable quantum phases. The first 
is a metastable Mott-insulating phase at (/>o = 0i = (t>-i = 0. 
The second one is a metastable transverse polar superfluid 
phase at 0i = 0_i — 0.38 and 0o = 0. 

Finally, we have also analyzed the ferromagnetic case. The 
presence of a magnetic field changes the energy functional of 
the ground state in such a way that first-order phase transitions 
and metastable phases are possible. 

IV. DISCUSSION AND CONCLUSIONS 

We have analyzed the ground-state phase diagram for spin- 
less and spin-1 atoms in period-2 superlattices. The dynam- 
ics within the unit cells was included exactly and the tunnel- 
ing between unit cells in a mean-field approximation. We 
discussed several methods to treat this mean-field Hamilto- 
nian and concluded that in the spinless case a simple stabil- 
ity analysis is sufficient to determine whether the system is 
Mott-insulating or superfluid. Using this method, we have 
first calculated the phase diagram for spinless bosons in opti- 
cal superlattices. In agreement with previous studies |22], we 
found a contraction of Mott lobes to loops for specific values 
of the energy offset. We have presented a detailed study of 
the various Mott phases which emerge when the chemical po- 
tential and the energy offset are varied. Furthermore, we have 
calculated the occupation numbers of single sites and found 
single-atom resonances in the Mott-insulating regime. These 
were known for isolated double-well potentials and were gen- 
eralized in this paper to the case of superlattices. We found 
clear fingerprints of single-atom resonances also in the den- 
sity of condensed bosons. 

In the case of spin-1 atoms the mean-field Hamiltonian 
shows a much richer quantum phase diagram. For anti- 
ferromagnetic interactions all Mott lobes are elongated to- 
wards higher tunneling amplitudes. Mott lobes with an even 
number of atoms at each lattice site are especially favored be- 
cause their atomic spins can couple to form spin singlets. A 
small, non-vanishing order parameter leads to increased atom 
number fluctuations and higher spin configurations, and, as a 
consequence, to higher on-site repulsion. Thus, the ground- 
state energy for small values of the order parameter is in- 
creased. For certain parameter regimes this leads to the ap- 



pearance of two local minima of the ground state energy func- 
tional (one at vanishing order parameter, one at a finite value 
of the order parameter) separated by an energy barrier. The 
higher one corresponds to a metastable quantum phase. Thus, 
the system shows a hysteretic behavior and the phase tran- 
sitions are of first order, whereas they are strictly of second 
order for the spinless case. For a realistic value of the spin- 
dependent interactions for sodium it depends on the parity of 
the atom number if the phase transition of a specific Mott lobe 
becomes first order or remains second order For smaller val- 
ues of the spin-dependent interactions all Mott lobes show 
first-order phase transitions, contrary to the case of spin-1 
atoms in usual lattices where the phase transition of the Mott 
lobe with one atom per site remains second order for all val- 
ues of the spin interaction fT?]. Because of the richer proper- 
ties of the energy functional in the spin-1 case it is no longer 
possible to determine the quantum phase of the system with 
a stability analysis of the Mott phase only. The ground-state 
energy functional for each point in parameter space has to be 
analyzed, and we have given a detailed analysis of the size of 
the various Mott phases in the atomic limit and pointed out 
the differences to the case of spinless bosons. We have also 
compared the single-atom resonances for spin-1 atoms with 
the case of spinless atoms and concluded that spin-dependent 
interactions change the occupation numbers of individual lat- 
tice sites. Spin-1 atoms in optical superlattices are therefore a 
model for mesoscopic magnetism. 

Finally, we have discussed the effects of magnetic fields 
by using an effective Hamiltonian which includes a quadratic 
Zeeman shift. For anti-ferromagnetic interactions magnetic 
fields break the degeneracy between different polar super- 
fluid phases. This leads to new classes of metastable phases 
and thus an even richer phase diagram. In the ferromagnetic 
case magnetic fields cause first-order phase transitions and 
metastable phases. These results apply to spin-1 atoms in su- 
perlattices as well as in usual lattices. 

In conclusion, we have shown that spinor bosons in optical 
superlattices show a rich phenomenology of different phases. 
We hope that our analysis will be a useful guideline for future 
experiments. 
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